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£ — | Small M -theories unify various models of a given family in the same way as the M-theory unifies a 

variety of superstring models. We consider this idea in application to the family of eigenvalue matrix 
models: their M-theory unifies various branches of Hermitean matrix model (including Dijkgraaf- 
Vafa partition functions) with Kontsevich r-function. Moreover, the corresponding duality relations 



o 



look like direct analogues of instanton and meron decompositions, familiar from Yang-Mills theory. 



2h 1 1 Introduction 



String theory 1 approach to analysis of QFT models involves the following seven steps. 



> 
• r- 1 

x. 

1. When one introduces particular QFT model it is usually defined as a perturbation expansion 
- ' around a given "classical" background. 

2. Perturbation theory is usually defined in terms of correlators that are Gaussian integrals. This 
leads to two corollaries. The nice one is the Wick theorem and the powerful diagram technique. The 
bad one is divergence of the perturbation series, which is asymptotic series due to the standard Dyson 
argument. This is causes by expanding integrands without respect to their asymptotics and, therefore, 
non-trivial phase/branch structure. 2 

3. The divergent perturbation series can be often presented in the form of a "genus expansion" 
asymptotic series in a single "Planck-constant" parameter, particular terms of the genus expansion 
being convergent multiple series in numerous other perturbation parameters (called coupling constants 
or time- variables). 

a E-mail: al@itep.ru 

b E-mail: mironov@itep.ru; mironov@lpi.ru 
c E-mail: morozov@itep.ru 

1 We use this term in the spirit of ref.Q, as a name for the broad modern version of quantum field theory (QFT) 
with numerous applications, of which various superstring models of sub-Planckian physics are important, but particular 
examples. 

Some more clever non-Gaussian expansions where weights are dictated by dominant non-quadratic terms in action 
are already in agenda (see |5] [3] for different conceptual approaches to the problem) , but are still far from any practical 
significance. 
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4. Non-perturbative partition function gets contributions from a variety of "classical" backgrounds, 
each contribution dressed by its own perturbative expansion. 

At this stage we get a representation of non-perturbative partition function - i.e. the only quantity 
which has an objective physical meaning - as a peculiar triple expansion: over "classical backgrounds", 
asymptotical quasiclassical (genus) expansion near each background, convergent series in powers of 
coupling constants for each term of the genus expansion. Actually the boundaries between the three 
levels of hierarchy are not well defined: changes of integration variables easily mix the levels. Worse 
than that, given this triple expansion (of which we can usually calculate only the first terms), it is 
rather difficult to understand the structure of the globally defined non-perturbative partition function 
(beyond any doubt, it is a nice function, just with sophisticated analytical structure, with numerous 
branchings and, often, essential singularities, which make perturbation series numerous, divergent and 
difficult to sew). The next three steps are targeted at curing this situation. 

5. Of the three levels the highest - sum over backgrounds - is understood the worst. It is often 
formulated in the language of "instantons" : the idea is that this sum has its own hierarchical structure 
and infinitely many terms in the sum are somehow expressed through multi-instanton configurations 
(that probably can be build from the single- instanton ones by acting a few operator generators). A 
nice realizations of this idea exist at the level of ordinary quantum mechanics see jlj |3] and in Ad QFT 
in the supersymmetric context (where perturbative contributions artificially cancel) 0. 

6. It is well known, however, that instanton calculus is hardly exhaustive: the proper generators, 
if any, would rather be made from instanton elementary constituents - merons [7J|S]. 3 Also the role 
of merons in the distribution of roles between the sum over backgrounds and genus expansions, i.e. 
between the two upper levels of triple expansion, remains under-investigated (as almost everything in 
meron physics). 

7. Different kinds of expansions - emerging when one starts from different QFT models at step 1 
- at non-perturbative level are unified in a global partition function. Since it unifies a priori different 
models it deserves a special name: nowadays it is often called M-theory. Being a true result of 
functional integration, the global partition function inherits a hidden symmetry ^01 (coming from the 
freedom to change integration variables and providing a non-trivial generalization of the Ward-like 
symmetries) : it belongs to a narrow class of integrable r-functions , satisfying physically- mysterious 
quadratic relations |1()|I12|. Various models, arising in different limits in the space of couplings (when 
one or another background gives dominant contribution), are said to be related by duality: sometime, 
but not always, duality transformations form a small group (but not obligatory Z<i or SL(2), some 
dualities can look like trialities etc, generically there is no small group at all). Switching on more and 
more perturbations, one enlarges M-theory, and at the end of the day one can arrive at the all-unifying 
"theory of everything" (nicknamed String Theory in |I]), where every two QFT models are "dual" 
in the sense that complicated correlators evaluated in one of them are converted by an appropriate 
analytical continuation into other complicated correlators in another one. However, one can stop 
at intermediate stages: by restricting classes of correlators (for example, to those of polynomials of 
original fields ), one obtains M-theories for sufficiently small varieties of QFT models. 

The task of the present paper will be a brief introduction into M-theory of ordinary (eigenvalue) 
matrix models This presentation is based on our recent paper ^3], which provided a unified 

derivation of significant results about relations between a priori different matrix models (Kontsevich 
model |14j , complex-matrix model |15| IT6] and Hermitean model |17( I18j with its many phases JH] > 
including the Dijkgraaf-Vafa |2(J| and other [2^ branches), obtained during the last decade |22|-|26| 

3 For Yang- Mills theory it is a puzzling problem: Yang-Mills instanton can be described as a process of creation and 
annihilation of a monopole-antimonopole pair, splitting of instanton into two merons is used to separate contributions 
of these individual quasiparticles and can provide proper variables to describe monopole condensation. Thus, the meron 
calculus should play the central role in confinement mechanism |7j. At the same time there is still no clear room for 
merons in the ADHM construction jjj]. 
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and based, in turn, on the long studies of eigenvalue models, since they were first introduced by 
F.Dyson in 1962 |T7] . 

The duality relations connect the partition functions of the Hermitean and Kontsevich matrix 
models. The third important eigenvalue model in the same M-theory is less popular but still important 
complex matrix model with partition function Zc(t). Actually, duality expresses Zq through Zk |13| . 

Like the Hermitean model, the complex one has many phases/branches, differing by shifts of the 
time- variables, like (13) and (15). We reserve the notation Zc{t) for the simplest, Gaussian phase 
(similarly, in the present paper Zk(t) denotes the simplest, Kontsevich branch of the GKM partition 
function; for dualities between different branches of GKM see |34| ) . 

Other models, especially the very important unitary one [27], can probably be also included into 
this matrix-model M-theory. 

2 The WEB of matrix-model dualities 

2.1 Basic instanton and meron decompositions 

Partition function of the Hermitean matrix model has many branches, some are labeled by two func- 
tions W(£) and /(£) [El EI]: Z\yj(t). These phases are distinguished by existence of a special - 
t'Hooft's - genus expansion, where genus is indeed the genus of the Feynman-t'Hooft fat-graph dia- 
grams. The leading (spherical) terms of this expansion are known as Dijkgraaf-Vafa partition functions 
[271] , The best studied is the Gaussian phase, Zq{€) = Zi M ^ 2 iMS (t), associated with W(£) = ^M£ 2 
and /(£) = 4MS = const. Here S = Ng, N being the size of the Hermitean matrix. Zoit) is the 
Toda-chain r-function ^H], while other Zwj(t) are r- functions of less known hierarchies. 

For polynomial W and /, Z\yj can be expressed through the Gaussian Zo(t), like the multi- 
instanton partition function would be expressed through the single-instanton one. If the derivative 
W(£) is a polynomial of degree n in £, W'(£) = rG=i(£ ~~ &)j 

Z WJ (t) = U wlG (t\tM,. . . *<»>){ ®U 4\t (i) )} (!) 

where time-variables t^' are obtained from t by shifts, so that they describe expansion around the 
extremum £j of W(£), and parameters Mj and iVj of Zq are made from W"(£j) and /(&), see |241 
Es3EHlE3[^ for the details. Operator . . . t^) intertwines between elementary constituents, 

such intertwiners will appear in all duality relations. All these operators have a very special form, 
typical for the Weyl-Moyal operators |321 13U| and their Batalin-Vilkovisky generalizations [33] : they 
are quadratic exponentials of the time-variables t and/or time-derivatives. 

However, the decomposition formula (^Q) is only the beginning. Like Yang-Mills instanton is made 
from two merons, the Gaussian partition function Zq{€) decomposes into the two Kontsevich r- 
functions Zk(t) [T5] : 

Z G (t) = U G \ K (t\T ± ){z K (T + ) Z K {r-)} (2) 

where the variables r± can be expressed through t in a variety of ways (for particular t— t± relations this 
formula was first suggested in ^2] and carefully derived in [22] )• Accordingly the meron decomposition 
exists for the other branches described by the polynomial W of degree n + 1: 

Zwj(t) = U wlK (t\T^\ . . . ,t£">){ C=i (4V?) ® 4V W ) } (3) 

4 To avoid a confusion, we remind that the term "M-theory" was introduced in I28| for a unified description of the web 
of superstring models. Progress with this superstring M-theory, revived interest |29| to the old landscape-style views on 
string theory (e.g. in |3UlH|) and finally brought more attention to the string theory per se, irrespective of its particular 
applications - even that great as unification of all fundamental interactions and of all cosmological scenarios. The M- 
theory of matrix models, considered in the present paper, has no direct relation to the superstring A/-theory (which also 
uses the matrix-models technique |31p. Though mathematically and conceptually similar (and this is what makes our 
study important), they concern absolutely different physical theories, the very small one, even with no time-evolution 
(but allowing a detailed and rigorous description), and enormously big one, big enough to include all our world along 
with entire variety of all other worlds, thinkable and unthinkable. 
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However, in variance with ©, this formula is not fully investigated in |T3], operator U\y\K is n °t yet 
explicitly constructed and the relation 



Uw\g®UU { g} k = Uw\k (4) 

also requires direct derivation. Entire theory of intertwining Weyl-Moyal operators and underlying 
♦-structures ^3] remains to be built. 

2.2 Polynomial representations (moment variables) 

Z G (t) and Zk(t) are themselves invariant under some restricted changes of variables |13j : 



Z G (t) = U G[G (t\t)Z G (t) (5) 



and 



Zk(t) = U k \ k (t\t)Z k (t) (6) 

In these formulas intertwiners U G \ G and Uk\k do not contain time-derivatives: they are just functions 
of time-variables and we do not write hats over them. As in all previous cases we do not write 
the somewhat lengthy formulas for transformations of time variables explicitly (they can be found 
in ref.^3]). I n © an d © these transformations are linear and depend on two free parameters. 
However, these few parameters can themselves be made arbitrary functions of times, and the two- 
parametric families of linear transformations generate infinitely large (still special) families of non- 
linear transformations. Moreover, among these non-linear transformations there is one that greatly 
simplifies Z G (t) and Zx{t): namely, all except the first two terms in genus expansions can be made 
polynomial in time- variables. We call these special choices i and f , and actually 



and 



Z G (t) = Z G (t) (7) 



Z k (t) = Z K (f) (8) 



In these formulas Z is defined in the following way: if logZ = ^2 p >og 2p ~ 2 F p is the genus expansion 
for Z, then log Z = X) P >2 9 2p ~ 2 F p and, all F p with p > 2 are polynomials of t or f variables (with 
powers growing as p increases). Since in © and © intertwiners are functions, and - as a corollary - 
in this case contain only terms with g~ 2 and g° in the exponent, they do not show up in © and ©. 
These formulas are the best possible illustration of statement 3 in the Introduction: contributions of 
given genus are not just convergent, in appropriate coordinates they are polynomial. 

The only word of caution to be made is that while © holds for conventional (t'Hooft's) genus 
expansion, in © another, - loop - expansion is implied, see [I5|. The polynomial coordinates exist the 
in Gaussian case for the t'Hooft expansion as well (3^1, associated variables are known as moments, 
but unlike © the corresponding relation 

Z G (t) = U G \ ACKM {t\n)Z ackm(h) (9) 

involves a new partition function Z ackm(iJ')- contributions of genera p > 2 to log Zackm{h) are 
polynomials in [i, but Zackm and Z G are different functions. 5 

To avoid possible confusion, we give a trivial example illustrating the difference between © and (j^J: 

■ - 1 / -l\ 
sin t — sin t — — {ft — ft ) 

where i = t + 2n and fi = e lt . The first equality is a model of with a single function Z(t) = sin(t), while the second 
equality relates two different functions, sin a; and ji{x — x^ 1 ). 
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3 Group theory approach 



Partition functions of eigenvalue matrix models are best defined as the highest weight vectors in the 
special Verma modules of Virasoro algebras, i.e. are eigenvectors of their maximal (Borel) subgroups: 

L l _Zi = (10) 

Here I = G, K,C and Virasoro loop operators are: 

00 dz 2 [ °° d d 2 

n=-l Z [k=0 ° lk+n a+b=n OTa ° Tb 

oo i j f oo o r>2 I | i 



Z'' 

n=-l I fc=0 



00 ,7,2 oo ^ * 



In these formulas 

4 = -~M5 k>2 + r£ = --M<J M + r fe (12) 

Zyj/j satisfies the equations 

L W Z WJ = (13) 
where L_(z) differs from L_(z) only by a shift: for generic W(£) it can be written as 

£ = - W«) + E fc^^- 1 (14) 

k k 

The function /(£) parameterizes different formal-series solutions to this loop equation, possessing 
t'Hooft's genus expansion. In the Gaussian case of W(£) = \M£ 2 , it reduces to /(£) = AMS = const. 

It is expected that similar representations exist for generic non-perturbative partition functions 
(though even the proper analogues of Virasoro algebras are not yet known in most cases). Thus, the 
relevant formalism of non-linear realizations of symmetries is of universal importance. The realizations 
are non-linear because of the shifts of time-variables in the Virasoro equations. To understand the 
structure and symmetries (dualities) of eigenvectors in non-linear representations, it deserves to begin 
from the simplest examples, with ordinary operators instead of the loop ones. We return to matrix 
models in section |SJ 

4 Examples of non-linear realizations of symmetries 

In this section x, y, z, t, r will play the role of the time-variables t and r of the previous sections. 
4.1 E* <g> SO (2) group 

Consider the SO(2) invariant functions Zjq on the real plane M 2 that realize the weight N represen- 
tation of the group of multiplications by real numbers, 

y'7r- x '7r) z " = ° 

Vj (15) 

x'^- + y'^]z N = NZ 
ox ay ) 
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where the shifted variables x' = Rcos<p = x + r, y' = Rsincf) = y parameterize the real plane. Since 
there is no central extension, the maximal subalgebra in this case can be taken to coincide with the 
entire algebra. 

Now one may use the two different strategies to find Z^. First of all, one can just note that the 
SO{2) invariant functions on the real plane M 2 depend only on the radial variable R, while the second 
equation in ()15j) . the homogeneity condition implies that Zjy ~ R N . The main subtlety here is related 
to the shift of coordinates (x',y') — > (x,y) so that one should re-calculate the linearly realized SO(2) 
symmetry in coordinates (x',y') (</> — > </> + a) to its non-linear realization in coordinates with shifted 
origin (x, y). 

Thus the "partition function" 



Z N = const • R N ~ ( 1 + — + ^X^- ) (16) 
\ r r z I 

The other possible way of calculating is to construct the solution to ()15|) as a power series in x 
and y, normalized to start with unity: 

, . / 2x x 2 + y 2 \ N/2 N N(N-l) , N 2 n „s 

Z N (x,y)=[l + - + —^^ =l + V x+ , x 2 + — 2 y 2 + ... (17) 

Coefficients in this series can be recurrently determined from (|15j) (this is the only way to construct 
solutions available in matrix models). 

To compare these two approaches, one can note that there is the identity: 

Va Z(x, y) = z(^ — r(l — cos a) + x cos a + y sin a, — r sin a — xsina + y cos a) (18) 

Let us check it with for the power series: 

Nf . \ N(N-1)( 
H 1 — r(l — cos a) + xcosa + y sin a J H k 1 — r(l — cosa)+ 



N/2 



r 



2r 2 



+x cos a + y sin a) + ^-^^ — rsina — xsina + y cos a) +...= (19) 
N N(N-l) n N 9 

where both sides should be understood as power series in x,y and a. For example, the terms without 
x and y are 

-N(l - cos a) + -N sin 2 a + }rN(N - 1)(1 - cos a) 2 + . . . = 

^ ^ ■ 2 , ■ 2 « , ■ 4 « \ ^V, • 4 a \ (2,,! 
= — sin 2 a - 4 sin 2 4 sin 4 — {-...) H 4 sm 4 — + ...+... 

2 V 2 2 / 2 V 2 



For a = f one gets 



Z N (x,y) = Z N (-r + y, -r - x) (21) 



This is trivially true, if one knows that Z is a function of (r + x) 2 + y 2 , but is a non-trivial identity 
for formal series, if one does not know their sums explicitly. 

4.2 T-duality 

Much closer to the case of matrix models is the following example 6 : 



3 A simple one-dimensional example that is very close to the matrix model case is considered in detail in I3(i| . 
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By a change of variable t = e tT it is transformed into the heat/Shrodinger equation (among other 
things, this t — r relation emphasizes "un-naturalness" of expansions in integer powers of t, the 
expansion around the point t = is not very nice from the point of view of the heat equation): 

fa a 2 \ 

{-'Tr + a7<) z = (23 > 

with the generic solution 

Z(x\t) = J e ipx e ip2T c(p)dp = J t p2 e ipx c(p)dp (24) 

In this 2-variable example it is obvious that (|22|) defines Z(x\t) ambiguously, and solutions are pa- 
rameterized by a function c(p): a counterpart of /(£) in Zy/j- 

Particular solutions of the heat equation include: 

• the heat kernel on a complex line = real plane 



Z x (x\t) = y -^e~ lx / 4t , with ci(p) = l; (25) 
the periodic Jacobi ^-function (actually, this 9 = 8qq, but we omit the index 00 in what follows) 

oo oo 

J2 e m 2 r e mx ) ^ ^ = g tf(p_„); (26) 





( X 




|r) = e 




I) 









the heat kernel on a torus 



„ oo 



— e~ iX ^""^ ^—i-K 2 n 2 /T^mnx/T (27) 



r / I t 

1 n=— oo 



IX 



2V- 



IT 



Function cs(p), associated with Z3, actually coincides with C2(p): this is the corollary of T-duality, 
which states that 

Z 2 (x\t) = Z 3 (x\t) (28) 

Usually T-duality is treated as a transcendental relation, derived, for example, with the help of Pois- 
son re-summation trick. However, one may represent all solutions to the heat equation, including Z\, 

Z2 and Z3, as linear combinations of Hermite rescaled polynomials He n (x|r) = {—it) 11 / 2 H n 

(where we denoted through H n the standard Hermite polynomials) 

Heo(x|r) = 1, Hei(x|r) = ix, He2(x|r) = 2ir — x 2 , He3(x|r) = — 6xr — ix 3 , 

He 4 (x|r) = -12r 2 - 12ix 2 r + x 4 , ... 

since they solve the heat equation. Then, the relations between different partition functions become 
involved but elementary algebraic relations. 
For instance, in representation (|24j) 

OO 72 

etfrjpx = J- P-He n (x|r) (29) 

n=0 

while, in order to get an expansion of the dual representation into the Hermite polynomials, one needs 
to shift the variable r — > r 1 = T + r and expand in inverse powers of the background T. First, we 
expand in this way Z\. 

— e- 2 ^ — — ( 1 + -) ~ 1/2 exp { -i ( ^ i ( 1 + ir 1 ] - T im He2m(X ' T) (31) 
This formula is a particular case of the celebrated Moller formula, eq.(10. 13.22)] 
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Now for the dual representation, one needs to shift in this formula x by q to obtain 8 



_J_ e ~ix 2 /4r e iqx/r e -iq 2 /r 



E 

m>0 



He 2m (x - 2q\r) 



Am'T"m+l/2 



ml 



E > 

k,m>0 



m—k 



(2q) k r(m+i) He2m _ fc ( 



XT) 



T m +V2 p^fj fc!(2m-fc)! 

(33) 



At the level of the heat equation, the T-duality states that 

/ e -ixV4r\ - 1 / g #2 \ e -ix 2 /4r 3 

I a/t I \ 8t dx 2 J y/r dr 



+ 



d 2 



IX 



a 



x=const 



dx 2 T dx 



1 



d 



+ 



d 2 



i ^ J z/r=const <9^x/V^ 

Accordingly, for every solution Z(x\t) of the heat equation, the T-dual function 



Z(x|r) 



-ix 2 /4:T _. 



(34) 



(35) 



is also a solution. 



5 Spectral-surface approach 

In the case of the Virasoro constraints only a non-linear realization of the symmetry algebra is known: 
identification of coordinates, where it is realized linearly, remains a puzzle. It is known as the problem 
of true variables in string theory, and it remains unresolved at all levels, from quantum gravity and 
entire string theory to matrix models. Thus dualities, described by eqs.Q-© at our present state of 
knowledge remain non-trivial - as would be 1)21(1 at the level of (|19jl. without (|16ft . Therefore, these 
dualities should be somehow derived. In ref.^S] we did it with the help of an additional structure: 
by considering various loop operators LL as different asymptotics of a single global loop operator 
defined on an auxiliary spectral Riemann surface. The same technique is believed to be adequate for 
all M-theories, one just needs to identify appropriate spectral manifolds... 

Definition of the global loop operator actually involves several substructures Three are of 

universal importance, but we formulate them in matrix-model terms, where the spectral parameter z 
is a single complex variable: 

- the choice of the spectral curve X, where the loop parameter z is taking values (z does not lie 
just on the complex plane, rather z £ T,, and S can be a non-trivial complex curve, a ramified covering 
of the bare curve So, which can be either C or Riemann sphere C or complex torus C*); 

- the choice of the Krichever-Novikov-type [HE] algebra, to which the loop operators L(z) belong, 
it is specified by the choice of the structure algebra (the Virasoro algebra in the case of eigenvalue 
matrix models) and allowed singularities on a bare spectral curve So (in the Virasoro case allowed are 
poles and square-root singularities, for W n algebras ramifications of orders that are divisors of n are 
allowed), conversion to S resolves the ramifications so that only poles are allowed in z £ S; 

- the choice of a Borel subalgebra L_(z), which can have non-trivial eigenvectors, in terms of the 
spectral manifolds this minus-projection can be represented by an integral operator, 

L-{z) = j> K(z,z')L(z'). 
8 To derive this formula one suffices to note that, for the Hermite polynomials, 



(36) 
(32) 
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Particular operators L_(z) arise from the global L-{z) in vicinities of particular singularities, poles 
are associated with L_, quadratic ramifications - with and L_ (more general GKM r- functions 
would arise if the structure algebra was W n with higher order ramification points). Duality relations 
(EQ)-© arise when the contour C in (|36|) is decomposed into combinations of contours, surrounding 
different singularities. Most analytical changes of the spectral variable z — > z(z), dz/dz = break the 
structure of minus-projection, but some - preserve it, and they give rise to identities like ©-©• 

In the case of eigenvalue matrix models there are two additional substructures ^3], which at least 
simplify the theory, but can also appear to have more general significance. They reflect the existence of 
free- field formalism on Riemann surfaces, and thus are relevant to non-perturbative partition functions 
(r-functions) associated with 2d conformal models. 

- The Virasoro algebra can be embedded into the universal enveloping of the Kac-Moody algebra 
by the Sugawara construction [22], and this embedding, L{z) =: J(z) 2 : can be continued to Krichever- 
Novikov-type algebras (in the matrix-model case, the simplest Kac-Moody algebra £7(1) is relevant, but 
everything should work for non-perturbative partition functions, associated with all other Kac-Moody 
algebras - this is the most straightforward generalization of our construction); 

- An additional Seiberg-Witten structure ^Hl^D; a 1-form Qsw( z ) with special properties, can 
be used to construct the kernel 

A = TT^TTsi 9 ^) 4>V)) = n S w(z)k(z, z') (37) 

"SW(2 J x ' 

in ()36j) from the free-field correlator on the spectral curve; actually this kernel can be used to define 
a multiplication operation of the Kac-Moody currents 

(j*j)(z) = £k(z,z') : J(z') 2 : (38) 

which plays a big role in "string field theory" [12] and/or in "background- independent" description 
|43| of the matrix model M-theory along the lines of refs. |44| 145 j. 
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